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Abstract 

We analyse the events with two same-flavour, opposite-sign leptons + E^^^^ 
+ (jets) as expected in pp collisions at LHC within the framework of the minimal 
Supergravity Model. The objective is the determination of the parameters mo and 
mi/2 of this model (for a given value of tan/3). The signature -\- EJp"^^^ + (jets) 

selects the leptonic decays of X2 ~^ Xi'^^ i X2 ~^ ^lr^^ ~^ Xi^^l - We exploit 
the fact that the invariant dilepton mass distribution has a pronounced structure 
with a sharp edge at the kinematical endpoint even in such an inclusive final state 
over a significant part of parameter space. We determine the domain of parameter 
space where the edge is expected to be visible. We show that a measurement of 
this edge already constrains the model parameters essentially to three lines in the 
{mQ,mi/2) parameter plane. We work out a strategy to discriminate between the 
three-body leptonic decays of X2 ^^'^ the decays into sleptons Il,r- This procedure 
may make it possible to get information on SUSY particle masses already with low 
luminosity, Lint = 10'^ pb~^. 



1 Introduction 



If 'low-cncrgy' supersymmctry (SUSY) is realised in Nature it should show up at the 
Large Hadron Collider (LHC). Strongly interacting particles as gluinos and squarks will 
be most likely the first SUSY particles to be seen at LHC. Gluinos of mass less than ~ 2 
TeV and squarks of mass less than ~ 1.5 TeV [1-3] can be detected, covering in such a way 
the whole theoretically motivated parameter space. LHC is also a good laboratory for 
the search of electro- weakly interacting particles, e.g. sleptons [4,5]. In a recent paper [6] 
it was shown within the minimal supergravity (mSUGRA) [7] model that sleptons in the 
mass range of ~ 100 to 400 GeV can be detected at LHC by investigating the signature 
two leptons + E^'^^'^ + no jets. However, this final state where direct production (Drell- 
Yan) of sleptons predominates requires high luminosity. Lint — 10^ pb~^. Sleptons can be 
also produced indirectly in the decays of charginos and neutralinos, especially in X2 ~^ Il,rI 
decays. The charginos and neutralinos, can in turn be produced directly or come from 
gluinos and/or squarks. This leads to final states with > 2 leptons + E^^'^^ + {jets). 
Actually, this indirect slepton production through g, q decays has the largest cross-section 
in a sizable region of the parameter space accessible at LHC and could allow sleptons 
to be already revealed at L^nt = 10^ pb~^, i-e. simultaneously with strongly interacting 
sparticles. Thus, indirect production of sleptons can be more important for observing a 
slepton signal than direct one [8]. Moreover, in such a way the mass reach for sleptons 
search can be extended up to mj^ ~ 740 GeV. 

Having evidence for SUSY at LHC, one of the next tasks will be to find out the un- 
derlying model and to determine the model parameters. In this paper we work out a 
method to determine SUSY parameters and suggest a strategy for getting information on 
masses of SUSY particles by means of the signature + + [jets). Our study 

is made within the framework of the minimal supergravity model (mSUGRA) [7]. In 
this model all scalar particles (sfermions and Higgs bosons) have a common mass mo 
at Mgut ~ 10^^ GeV. The gaugino masses Mi,M2,M3 (corresponding to U(l), SU(2), 
and SU(3), respectively) unify to a common gaugino mass mi/2, and all trilinear cou- 
pling parameters Aijk have the same value Aq at Mgut- One also has unification of the 
electro weak and strong coupling parameters ai,i = 1,2,3 [9]. A further reduction of the 
parameters is given by invoking 'radiative symmetry breaking'. As a consequence, one 
has only the following input parameters: mo, mi/2, Aq, tan/3, sign(/x). Here tan/3 = ^, 
the ratio of the two vacuum expectation values of the two Higgs doublets, and fi is the 
Higgsino mass parameter. The whole SUSY particle spectrum can then be calculated by 
making use of renormalization group equations (RGE). This model is also incorporated 
in the Monte-Carlo generator ISASUSY [10] which is used in our analysis. 

This paper is aimed at determining the parameters mo and mi/2 with fixed tanp. 
Knowing these parameters, we can calculate the masses of the superpartners using RGE. 
For this purpose, we study in detail the leptonic decays of X2 which have some useful 
features. Within the mSUGRA model, X2 has two-body decays, X2 ~^ ^t,BF^ 1 the region 
mo ^ 0.5 • mi/2 of the parameter space, whereas in the region mo ^ 0.5-mi/2, mi/2 ^ 200 
GeV the x^ has three-body decays, ~^ l'^^~Xi- both regions the invariant dilepton 
mass spectrum (Mi+i-) has a maximum MJ+f!^, and therefore a pronounced structure with 
a sharp edge can be seen at the kinematical endpoint. This property was discussed first 
in ref. [11] in the case of three-body decays of X2 then in ref. [12] in the case of 
two-body decays. The generality of this feature, i.e. the observability of an edge in the 
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Mi+i- spectrum even in inclusive l^r^l'^~ and /+/~ + final states in a large part 

of the parameter space was shown in [13]. We will show how much the parameters are 
constrained by a measurement of the M'^'^^ value of the dilepton mass spectra. Moreover, 
we will discuss a method, based on the analysis of the Af;+/- spectrum, to find out whether 
the observed edge is due to the two-body or three-body decays of 

2 Sparticle masses in mSUGRA 

Within the Minimal Supersymmetric Standard Model (MSSM) the masses of the neu- 
tralinos are determined by the parameters M 'rni/2{Mz)^, yU, and tan(3 using Mi ~ 
^tan^Ow ~ 0.5M (Mi being the U(l) gaugino mass). In the following, we fix tanP — 2 
(we assume that tanj3 could be known from previous experiments) and take Aq — 0. 
In mSUGRA | fi \ quite generally turns out to be | /i |> M, so that in this case 
m(x2) — ""^(xf ) — 2 ■ 'rn(xi) ~ M. Both Xi X2 gaugino like, Xi is almost a 
pure B — ino, and X2 almost a pure W'^ — ino. In supergravity the slepton masses are 
given by [14]: 

= mo + 0.15mi/2 - sin'^ewMlcos2(3 (1) 
mf^ = mo + 0.52mi/2 - 1/2(1 - 2sin^ew)Mlcos2p (2) 

ml^ml + 0.52m?/2 + 1 /2M|cos2/3 (3) 

Analogous equations exist for squarks [14]. Therefore, when the parameters mo, mi/2 
and tan f3 are known we can calculate all sparticle masses. A special case is the third 
generation of squarks and sleptons, where L— R mixing plays a crucial role. 

3 Production and leptonic decay of X2 

Neutralinos X2 can be produced at the LHC through a Drell-Yan mechanism (direct 
production), in association with strongly interacting sparticles, or in the decay chain of 
gluinos and squarks (indirect production). Gluino and squark pair production processes 
are the dominant source of x^'s because of large strong interaction cross-sections. The 
branching ratios of gluino and squark decays into X2 are also sizable and are shown in 
fig.l as a function of the model parameters mo and mi/2. One can see that the regions 
in (mo,mi/2) plane where the decays g ^ X2 + X and g — > X2 + ^ are open, are comple- 
mentary. In the region mo ^ 1.47 ■ mi/2 gluinos are lighter than squarks and can decay 
into while squarks prefer to decay into gluinos. In the region mo < 1.47 ■mi/2 squarks 
are lighter than gluinos and can decay into x^^ but gluinos decay into squarks. Hence the 
decays g x^ + X (q ^ X2 + ^) ^nd q ^ g ^ X2 + ^ (9 ^ Q ^ X2 + ^) can coexist (see 
also fig.l). Fig.2(a,b) shows a x Br for indirect X2 production from gluinos and squarks 
as a function of mo and mi/2- 

There are three different leptonic decays of X2 interesting for our study: X2 ^^^~Xij 
X2 ~^ ^R^^ and X2 ~^ ^L^^- When decays of X2 to sleptons are allowed, sleptons decay 
directly into the lightest supersymmetric particle (lsp= Xi) with Br(/J^^ Xi^^) — 100%. 
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Branching ratios of X2 decays into leptons directly and via sleptons are shown in fig. 3 as 
a function of mo and mi/2. The regions where these decays are kinematically allowed 
are complementary in the parameter plane, depending on whether X2 is lighter or heavier 
than Il,r- Thus, one can distinguish three domains in the (mo,mi/2) plane, which are 
(also see fig.4): 

domain I (mo ^ 0.5 • mi/2, 1^1/2 ^ 200 GeV ) : m^o < m^-^ ^ , with X2 ^ Xi ^ ^ 
domain II (0.45 • mi/2 ^ mo < 0.5 • mi/2 ) : < < mj^ , with X2 ^ 

domain III (mo ^ 0.45 • mi/2 ) : m^o > mj^ , with X2 ^ k^ 



In domain III, the decay X2 ~^ would also be kinematically allowed, but since the 
B — ino component of X2 is very small, the coupling to IrI is also small. Therefore, the 
decay X2 ~^ is very much suppressed in the whole domain. 

In fig. 5 we show the regions for a x Br{x2 Xi^~^^~) ~ 1 ^iid 0.1 pb, a x Br{x2 ~^ 
IrI"^ X\l^l~) - 1, 0.1 pb and a x Br{xl ^ ^ X?^^^") 1, 0.002 pb in the 
(mo, mi/2) plane from indirect and associated x^ production followed by decays to x?/^/~ 
final states directly or via sleptons. One can see that there are regions in domains II and 
III where the mentioned decays coexist. Finally, in fig. 6 we show a x Br for indirect and 
associated production of x^ decaying into leptons directly or via sleptons. 



4 Determination of mo, mi/2 and sparticle masses 

In order to determine the parameters mo and mi/2 we will exploit in the following the 
kinematical features of the two- and three-body leptonic decays of x^. 

As pointed out in [11], the decay X2 ~^ ^~^^~Xi (domain I) has the useful kinematical 
property that the invariant mass of the two leptons Mi+i- has a maximum at 

M,'?- = m^o-m^o, (4) 

whereas for the decays X2 ~^ ^l,r^^ ~^ l^l~Xi (domain II and III) the maximum of 
Mi+i- is given by [12]: 

./(r»^„ - niDinif - m|o) 

j^max ^ V 

Thus, the M;+/- distribution has a very characteristic shape with a sharp edge at the 
kinematical endpoint M"^^-. 

As the main source of x^'s is their indirect production in gluino and squark decays, 
the most suitable signature for selecting the X2 decays is provided by the topology with 
two same-flavour opposite-sign leptons accompanied by large missing transverse energy 
and usually accompanied by a high multiplicity of jets. In this paper we thus concentrate 
on the two same- flavour, opposite-sign leptons + E'lp'^^^ -\- {jets) channel, where the final 
state leptons are electrons and muons. 
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4.1 CMS detector simulation 



The simulations are done at the particle level, with parametrised detector responses based 
on detailed detector simulations. These parametrisations are adequate for the level of 
detector properties we want to investigate, and are the only practical ones in view of the 
multiplicity and complexity of the final state signal and background channels investigated. 
The essential ingredients for the investigation of SUSY channels are the response to jets, 
i?™""*, the lepton identification and isolation capabilities of the detector, and the capability 
to tag b-jets. 

The CMS detector simulation program CMS JET 3.2 [15] is used. It incorporates 
the full electro-magnetic (ECAL) and hadronic (HCAL) calorimeter granularity, and in- 
cludes main calorimeter system cracks in rapidity and azimuth. The energy resolutions 
for muons, electrons (photons), hadrons and jets are parametrised. Transverse and longi- 
tudinal shower profiles are also included through appropriate parametrisations. The main 
detector features incorporated in the Monte-Carlo description are: 

• Hadronic tracks, muons and electrons are measured up to | 77 |=2.4 

• Deflection of charged particles due to the 4 T magnetic fleld is included. 

• The resolution for the muon system is parametrised according to [16]. 

• The calorimetric coverage goes up to at | 77 | =5 for the HCAL and | r] |=2.6 for the 
ECAL. 

• ECAL energy resolution parametrized as: 

AE/E = 5%/^/E © 0.5% 

• HCAL energy resolution is parametrised according to [17] as a function of ry; a typical 
hadron resolution is: 

AE/E^80%/Ve®7% 

• Energy resolution for very forward calorimeter ( VFCAL) , in the parallel plate cham- 
bers option: 

AE/E = 90%/y/E © 3% 

• Granularity of calorimeters: 



V - 


range 


A(p X At] 


ECAL (barrel) 


1 V 1 
1.57 <| 


< 1.57 
rj \< 1.65 


0.015x0.015 
crack 


ECAL(endcap) 


1.65 <h |< 2 
2. <| ri \< 2.35 
2.35 <| T] \< 2.61 


0.022x0.022 

0.029x0.029 
0.043x0.043 


HCAL 


1 V 1 
2.26 < 
2.6 < 


< 2.26 
T] \< 2.6 

h l< 3 


0.087x0.087 
0.174x0.175 
0.195x0.349 


VFCAL 


3 <| 

4 < 


V 1 

V 1 


<4 
< 5 


10x10 cm^ 
5x5 cmP' 



• Ethreshoid on cells: 

- ECAL: Ethreshoid =50 MeV 
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- HCAL: Ethreshoid =250 MeV 

- VFCAL: Ethre shold 

=500 MeV 

• a modified UAl- jet finding algorithm with a cone size of AR=0.9 (for description 
see CMSJET 3.5 [15]) is used for jet reconstruction. 



4.2 Observability of edges in invariant dilepton mass distribu- 
tions 

In this chapter we determine the regions in the (mo, mi/2) parameter plane, where the 
characteristic edge in the Mi+i- distribution can be observed in inclusive final states with 
two same-flavour, opposite-sign leptons + + {jets) with different luminosities at 

LHC. 

Standard Model background processes are generated with PYTHIA 5.7 [18]. We use 
CTEQ2L structure functions. The largest background is due to tt production, with both 
1^'s decaying into leptons, or one of the leptons from a W decay and the other from 
the 6-decay of the same t-quark. We also considered other SM backgrounds: W + jets, 
WW, WZ, hb and rr-pair production, with decays into electrons and muons. Chargino 
pair production xfxi is the largest SUSY background but gives a small contribution 
compared to the signal. 

To observe an edge in the Mi+i- distributions with the statistics provided by an inte- 
grated luminosity L^nt — lQ^ph~^ in a significant part of the (mo,mi/2) parameter plane, 
it is enough to require two hard isolated leptons {p^'^ > 15 GeV) accompanied by large 
missing energy, i?™*'' > 100 GeV. Our criterion for observing an edge in the Mi+i- distri- 
bution contains two requirements: {Nev — Nb)/ ^/Nev ^ 5 and {Nev — Nb)/Nb ^ 1.3, 
where Nev is the number of events with Mi+i- < MJ+p, and Nb is number of the ex- 
pected background events. Fig. 7 shows the invariant mass spectra of the two leptons at 
various (7710,7711/2) points from domains 1, II and III, respectively. The observability of 
the "edge" varies from 77(7 and signal to background ratio 31 at point (200,160) to 27a 
and a signal to background ratio 2.3 at point (60,230). The appearance of the edges in 
the distributions is sufficiently pronounced already with L^^t — 10^ pb"^ in a significant 
part of (7710,7711/2) parameter plane, see fig. 9. The edge position can be measured with a 
precision of ~ 0.5 GeV. 

With increasing rriQ and 7/21/2 cross-sections are decreasing, therefore higher luminosity 
and harder cuts are needed. To achieve maximal reach in 7771/2 with L^nt = 10*pb~^ 
for points from domain III, a cut up to > 300 GeV is necessary to suppress the 

background sufficiently. For points with large 777o (domain I) the transverse momentum pt 
of the leptons and i?™ are not very large, but there are more hard jets due to gluino and 
squark decays. Thus for these points we keep the same cuts for leptons and missing energy 
as before {p^'^ > 15 GeV, > 100 GeV) and require in addition a jet multiplicity 

Njet > 3, with energy E^^ > 100 GeV, in the rapidity range | rjjet \< 3.5. To optimise 
the edge visibility we also apply an azimuthal angle cut, < 120°. For points 

from domain II, the jet multiplicity requirement is also helpful. Right sleptons are too 
light to provide large lepton pt and i?™*'', and to use cuts on pip and E^^'"^^ alone is not 
very advantageous. With L^^t = 10^ pb~^, to suppress the background at larger accessible 
7770 ) 7r7i/2 values, we have to require at least 2 or 3 jets, depending on the mo, mi/2 region 
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to be explored. Fig.8 shows invariant dilepton mass distributions at some {mo,mi/2) 
points close to maximum reach with L^^t — lO^pfe"^ and L^^t — 10^ pb"^ respectively. 

The regions of the (mo,mi/2) parameter plane where an edge in the Mi+i- spectra 
can be observed at different luminosities are shown in fig.9. In fig. 10 we show separately 
the three domains where an edge due to X2 ~^ ^^Xii ^R^ ^^'^ ^lI decays can be observed 
at Lint = 10^p6~^. One can notice a small overlapping region, where we expect to 
observe two edges, due to X2 ~^ ^^^'Xi to X2 ~^ ^r^^ ~^ ^^^~Xi decays {case 1). 
With increasing luminosity and correspondingly higher statistics, this overlapping region 
increases, see figs. 11 and 12. These plots show the same as fig. 10, but for Lint = lO'^p^"^ 
and Lint = 10^ pb~^, respectively. An additional region appears where two edges can be 
observed simultaneously, due to — >■ IrI^ — >■ l^l~Xi and it^^ ~^ l^l~Xi decays 

{case 2). These regions {case 1 and 2) are due to the coexistence of different decay 
modes has been seen in fig. 5. An example of a Mi+i- distribution for case 1 is shown in 
fig.l3. 

Therefore, to a given integrated luminosity at LHC {Lint = 10^ pb~^ to 10^ pb~^) 
there corresponds a definite parameter region where the characteristic structure in the 
Mi+i- distribution can be seen. This fact already gives a preliminary information about 
the parameters mo and mi/2. The observation of two edges would give even stronger 
constraints. 

4.3 MJ+f- analysis of (mo,mi/2) parameter plane 

Within mSUGRA all sparticle masses for every point (mo, mi/2) can be calculated. The 
expected edge position M™"5 in the dilepton mass distribution can then be obtained 
from eqs.(4) and (5). In fig. 14 we show contours for various expected values of MJ+f- in 
the (mo,mi/2) parameter plane. Different lines with the same value of belong to 

domains I, II and III (with the corresponding decay mode of x^). The region of MJ+f- 
becoming accessible at LHC is: 

for X2 ^ Xil^l~ - 50 GeV < MJ^I^ < 90 GeV 

X^ ^ IrI - M^-- > 10 GeV 

Xl IlI - M^p > 20 GeV 

More specifically, at low luminosity. Lint = 10^ pb at the beginning of the LHC 
operation, the accessible values of M/^"5 lie in the following ranges (see figs. 10 and 14): 

for xl Xil^l' - 50 GeV < MJI^f^ < 80 GeV 

X2 IrI - 10 GeV < MJl^p < 110 GeV 

IlI - 20 GeV < M^^^^ < 120 GeV 

It follows from the discussion above that a measurement of in the dilepton mass 

distribution, with a single edge, constrains the parameters in general to three lines in the 
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(mo,mi/2) parameter plane. In case of MJ^p > 90 GeV the constraint is stronger, there 
are just two possible hnes. The most favourable case is when the measured MJ^p value 
is large, MJ+'}^ > 180 GeV. Then one is left with a single hne in the (mo, m-1/2) parameter 
plane. For the present study we have chosen, as an example of the general situation, the 
case of MJ+p = 74 ± 1 GcV, with three lines corresponding to the domains I, II and III, 
respectively. The next step is to find out which line in the (mo,mi/2) plane is the right 
one. To this purpose we have analysed points along these lines, given in tables 1-3. 
The study is made for the low luminosity case. Lint — 10^ pb~^. 



Table 1. MJ^f^ values (in GeV) at the investigated (mo,mi/2) points from domain I, 

X^^X? + ^+ + ^- 



(120,160) 


(130,160) 


(180,160) 


(200,160) (220,160) (240,160) (290,160) (350,160) 


j^rnax 74 


74 


74 


74 74 74 74 73.7 



Table 2. MJ!^p values (in GeV) at the investigated (mo,mi/2) points from domain II, 



(mo,mi/2) 


(80,162) 


(90,170) 


(105,180) 


(110,187) 


(120,195) 


A/fmax 


74 


74 


73 


75 


73 



Table 3. M/^*"5 values (in GeV) at the investigated (mo,mi/2) points from domain III, 



(mo,mi/2) 


(20,195) 


(40,210) 


(60,230) 


(80,255) 




73 


73 


73 


73 



One should first notice that the observation of two edges at Lmt = 10^ pb~^ would 
determine the (mo,mi/2) point uniquely. This is due to the fact that the set of the edge 
position values in the Mi+i- spectrum is different at each point of the parameter region 
{case 1), where two edges are expected to be observed at 10^ pb"^, see figs. 10 and 14. With 
a luminosity Lint — 10^ pb~^ the positions of the two edges will fix two (mo, mi/2) points, 
belonging to domains II and III and corresponding to case 1 and case 2, respectively. 
At high luminosity. Lint = lO'"* pb^^, the observation of two edges can give up to three 
possible (mo,mi/2) points. One of them is from domain II, corresponding to case 1. The 
lines M;™"5 = const corresponding to X2 ~^ ^R^ decays have the form of an ellipse and can 
cross the MJ+J^ ~ const hnes corresponding to X2 ^ ^lI decays twice. Hence, two points 
with the same set of edge positions in the Mi+i- spectrum can be found in domain III 
corresponding to case 2. A discrimination between these points is possible on basis of the 
event kinematics, and/or by an analysis of the total event rate and the relative number 
of events corresponding to the two peaks, (see figs. 5, 10-12). 
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4.4 Discrimination between different X2 leptonic decays 

For a definite value of the edge position MJ+p one expects a different shape of the Mi+i~ 
distributions in two- and three-body decays (see fig. 13, where the first peak is due to a 
two-body decay of X2 the second one due to a three-body decay) . As we have seen 
from figs. 7 and 8 the signal events contribute in the interval < Mi+i- < MJ^f^. In the 
following wc only consider events in this mass region. The average value <Mi+i-> of 
signal and background events of this mass region is shown in fig. 15 as a function of niQ. 
The errors are calculated by taking into account the statistical error, a systematic error 
in the measurement of the edge position, and a systematic error of 30 % for background 
uncertainty (the main background is tt). One clearly sees that <Mi+i-> is significantly 
smaller in the case of a direct three-body decay X2 ~^ ^^Xi- Thus the shape of the dilepton 
mass spectrum already allows one to decide whether X2 decays into a slepton or not. 

In order to distinguish between domains II and III, we suggest to use the fact that in 
general the contour fines with the same MJ^i^ for right and left sleptons have no overlap 
in m 1/2 in regions of parameter space which are accessible at a given luminosity, see figs. 9 
and 14 as an example. It means that the masses of Xi's are different for these two lines, 
and hence i?™"'* is expected to be different. In fig. 16 we show the <i?™**^> values for 
events with X2 ~^ ^t,R^^ ~^ ^'^I'Xi decays after the cuts p^'^ > 15 GeV and E™*'''' > 100 
GeV, Mi+i- < M;+^5. The errors are calculated by taking into account the statistical 
error and a systematic error in the measurement of the edge position. As can be seen 
from fig. 16, is larger in the case of X2 ~^ ^l^^ ~^ ^^^~Xi than in the case of 

X2 ~^ ^R^^ ~^ l~^l~Xi ^ expected. 

4.4.1 Event rate analysis 

When the correct MJ+p line is chosen, the last step is to find the point {rriQ, 7711/2) on 
this line. In general the cross section falls with increasing mi/2 and mo. Thus, we study 
the event rate along the corresponding M^+"5 line. We first discuss the domain III, where 
the situation is simpler. For the event rate analysis at Lint — 10^ pb~^, to reduce the 
uncertainties due to background, wc use a harder cut on E'™''*, E"™** > 130 GeV. The 
dependence of the expected event rate on mo is shown in fig. 17. The errors are calculated 
by taking into account the statistical error and a systematic error of 30 % for background 
uncertainty. A systematic error due to tlic precision of the edge position measurement is 
also taken into account. From the observed event rate we can then determine mo with a 
good accuracy, Sttiq ~ 4 GeV. The parameter mi/2 is then given by the Af line in the 
(mo,mi/2) plane. The precision obtained in such a way is 61711/2 — 4 GeV. 

In domain II, the event rate along a line of definite Af™"5 is first increasing and 
then decreasing with mo, see fig. 18. Tliis is mainly due to the change in the branching 
ratios (see fig. 4). The dependence of the event rate is, however, such that mo cannot be 
determined uniquely. To a given event rate there correspond in general two mo values. 
The ambiguity can, however, be solved at high luminosity L^^t — 10^ pb~^, when two 
edges in the Mi+i~ distribution can be observed. 

For domain I, the mo dependence of the event rate is shown in fig. 19a, again for 
^max ^ 74 ± 1 GeV. Notice the steep increase of the rate at mo ~ 120 - 130 GeV. This 
is due to the fact that the decay channel X2 ~^ ^^Xi is just opening in this region. As can 
be seen from the curve in fig. 19a, there is an ambiguity in the determination of mo if the 
event rate is in the region 3700 < Nev ^ 5600 or 120 GeV < mo < 240 GeV. Here it 
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helps if we look at the average number of jets <Njet> in the events under study. Fig. 19b 
shows <Njet> as a function of mo- <Njet> is rising with mo as more jets are produced as 
the squarks become heavier. With the measured <Njet> we can resolve the ambiguity in 
the mentioned region 120 GeV ^ mo < 240 GeV and thus determine mo with 5mo ~ 7 — 3 
GeV. 

5 Conclusions 

In this paper we have performed a detailed analysis of events with the signature + 
^miss _|_ (^jets) to be expected in pp coUisions at LHC. Our aim has been to determine 
the parameters mo and mi/2 of the Minimal Supergravity Model and to get information 
on the mass spectrum of SUSY particles, assuming knowledge of tan(3 from previous 
experiments. Wc have exploited the property of the X2 leptonic decays X2 ~^ ^~^^~XiJ 
X2 ~^ Il,rI Xi that the invariant mass of the two final leptons has a maximum, 

clearly visible even in inclusive production. We have determined for different 
luminosities the regions in the (mo,mi/2) parameter plane where one or two edges can 
be observed in the invariant dilepton mass distributions. These regions already give 
preliminary information about the model parameters. The appearance of the edges in 
the Mi+i~ distributions can be already seen with a luminosity L^nt = 10'^ pb^^. Therefore 
we have concentrated on a low luminosity study. On the other hand, in case no such 
observation will be made at this luminosity, the corresponding parameter region can be 
excluded, and the same analysis can be done at higher luminosity. 

We have shown that a measurement of the M™°5 value constrains the parameters 
mainly to three lines in the (mo,mi/2) parameter plane. The lines correspond to the 

decay modes X2 ~^ l^l~Xi: ^lI — l'^l~Xi: ^rI I'^l^Xi respectively. We have worked 
out a method to discriminate the three-body from the two-body X2 decays. In the case 
of three-body X2 decays the parameter mi/2 can be determined by the measured value of 
MJ+p with a precision of ~ 0.5 GeV. The parameter rriQ can then be determined from 
the observed event rate with a precision of 7-3 GeV. In the case of two-body X2 decays, 
a measurement of the missing transverse energy can allow one to distinguish between the 
two possible decays X2 ~^ and x^ — > IrI, but a more detailed study is needed. By an 
event rate analysis along the corresponding line in the (mo, mi/2) plane we can determine 
mo and mi/2, ^trq ~ 5mi/2 ~ 4 GeV. 

Knowing mo, mi/2 and tan/3, the masses of all SUSY particles (except for the 3rd 
generation of squarks and sleptons) are calculated by RGE. The precisions which can be 
achieved are ~ 1 — 6 GeV. 

In such a way it is possible to obtain information about SUSY particle masses already 
with low luminosity {L = 10^ pb~^) even without having direct experimental evidence 
for their existence. This is especially important for sleptons in a parameter region where 
high luminosity would be necessary to detect them through direct production. 

This study has been performed for tanP — 2, but it is also possible for high values 
of tanp. Most likely, for large tanP > 30 a higher luminosity will be needed because of 
smaller branching ratios of the X2 Icptonic decays. 

Let us mention some further interesting aspects of this work. Selecting the two- 
body X2 leptonic decays by our method represents an indirect evidence for sleptons in 
the framework of mSUGRA. In this way it is possible to probe slepton masses up to 
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~ 740 GeV well beyond what is possible in direct [5,6] searches. As it has been shown 
in this study, the edge in the invariant dilepton mass distributions is expected to appear 
at M™"5 > 10 GeV, being quite generally a signal for a two- or three-body decay of 
some abundantly produced heavy object. Hence such an observation may serve as a first 
evidence for physics Beyond the Standard Model, and if observed with significant i?™'^'' 
it would be a clear evidence for SUSY, more specifically for X2 production. 
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Figure 1: Decay branching ratios as a function of mo and mi/2 (in GcV) for: a) g ^ X2+^) 
b) Ml ^ X2 + ^ 9 ^ + ^ , ul ^ 9 + ^ , for tan /3 = 2, Aq = 0, /i < 0. 
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Figure 2: Sigma times branching ratios as a function of mo and mi/2 (in GeV) for indirect 
X2 production from gluinos (a) and squarks (b), for tan (3 — 2, Aq — 0, /i < 0. 
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Figure 3: Branching ratios of X2 decays: a) X2 ~^ Xi^^^ ^ ^) X2 ^ ^l^^ c) X2 ~^ ^r^^ 
as a function of mo and mi/2, for tan P — 2, Aq — 0, /i < 0. 
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Figure 4: Domains of the decays X2 ~^ Xil^l~ (dashed hne), X2 ~^ ^t^^ (sohd hne) and 
X2 ~^ ^R^^ (dashed-dottcd line) in the (mo, "^1/2) plane, corresponding to decay branching 
ratios in excess of 1% and 10% respectively, tan/3 — 2, Aq = 0, /i < 0. 
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Figure 5: Contour lines for cross-section times branching ratios in the (mo, mi/2) plane 
for indirect and associated X2 production followed by decays: X2 ~^ xV^^^" (dashed 
line), xl II^^ ^ x\l^l~ (solid line) and X2 ^ ^ Xi^^l' (dashed-dotted line), 
tan/3 = 2, Ao = 0, /X < 0. 
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Figure 7: Invariant mass distribution of two same-flavour, opposite-sign leptons at various 
(mo,mi/2) points from domain I, II and III for L^^t — 10^ pb~^. SM background is also 
shown (dashed hne). 
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Figure 8: Invariant mass distribution of two same-flavour, opposite-sign leptons at 

(mo,mi/2) points from domain I, II and III close to the experimental reach at corre- 
sponding luminosities Lint = 10^ pb~^ and 10^ pb~^. SM background is also shown 
(dashed line). 
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Figure 9: Observability of edges in invariant dilepton mass distribution with luminosities 
10^ pb^^ (dashed line), 10'^ pb^^ (solid line) and 10^ pb~^ (dashed-dotted line). Also shown 
are the explorable domain in sparticle searches at LEP2 (300 pb~^) and the Tevatron (1 
fb~^), theoretically and experimentally excluded regions [19]. 
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Figure 10: Domains where the observed edge in the Mi+i- distribution is due to the 
decays ~^ ^l^^ ~^ Xil~^l~ (sohd hue), X2 ~^ ^r^^ ~^ Xi^^l^ (dashed-dotted hne), 
X2 ^ Xll'^l' (dashed hne), L^^t = 10^ pb-^ 
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Figure 11: Domains where the observed edge in the Mi+i- distribution is due to the 
decays ~^ ^l^^ ~^ Xil~^l~ (solid line), X2 ~^ ^r^^ ~^ Xi^^l^ (dashed-dotted line), 
X2 ^ Xll'^l' (dashed line), L^^t = 10^ pb-^ 
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Figure 12: Domains where the observed edge in the Mi+i- distribution is due to the 
decays ~^ ^l^^ ~^ Xil~^l~ (sohd hne), X2 ~^ ^r^^ ~^ Xi^^l^ (dashed-dotted hne), 
X2 ^ Xll'^l' (dashed hne), L^^t = 10^ pb-^ 
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Figure 13: Invariant mass distribution of two same-flavour, opposite-sign leptons at the 
point mo = 125 GeV, mi/2 = 170 GeV, L^^t = 10^ pb~^ SM background is also shown 
(dashed line). The edge at Mi+i- = 38 GeV is due to the two-body xl ^ IrI"^ ^ Xil^l~ 
decays and the edge at Mi+i- — 80 GeV is due to the three-body X2 ~^ Xi^"*"^~ decay. 
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Figure 14: Contour lines of expected MJ+f- values (in GeV) in the invariant dilepton 
mass distribution corresponding to the three different X2 decay modes in the region of the 
(777.0,777.1/2) parameter plane accessible with 10^ pb~^. 
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Figure 15: Average value of invariant dilepton mass, Mi+i- < M™^5 at (mo,mi/2) points 
with MJI^p = 74± 1 GeV corresponding to the decays ^ Xil^l', xl ^ ^ Xil^l~ 
and X2 ^ ^ X?^+^" as a function of mo; Pr" > 15 GeV, E^''' > 100 GeV, Li„t = 10^ 
pb-^ 
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Figure 16: Average value of at (mo, mi/2) points with MJ^p = 74 ± 1 GeV corre- 

sponding to the decays X2 ~^ ^l^^ ~^ Xi^~^^~ cind X2 ^ ^■bJ'^ ~^ as a function of 

tni/2, Lint = 10^ pb~^ Event selection criteria are: p^'^ > 15 GeV, > 100 GeV and 
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Figure 17: Expected 2l + EJp^'' event rate with L^^^ = 10^ pb"^ along M;^^5 = 74±1 GeV 
contour line in domain III (x" — ^ if^"^ ~^ Xi^^^~) ^ ^ function of mo- Event selection 
criteria are: pj'' > 15 GeV, E^'^^' > 130 GeV and Mi+i- < MJ^^ at corresponding 
points. 
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Figure 18: Expected 21 + Eip^'' event rate with = 10^ ph'^ along MJ^f^ = 74 ± 1 
GeV contour line in domain II Xi^^^ ) ^ function of itlq. Event selection 

criteria are: p^'^ > 15 GeV, E^*''* > 130 GeV and Mi+i- < MJ^p at corresponding 
points. 
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Figure 19: a) Expected 21 + EJp'^' event rate with 1^ = 10^ pb"^ along MJ^p = 74 ± 1 
GeV contour line in domain I Xi^~^^~) ^ function of ttiq. Event selection criteria 
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are: p^" > 15 GeV, E^^^'^ > 130 GeV and Mi+i- < MJ+p at corresponding points; 
b) Average number of jets (Ej^^ > 30 GeV, | rjjet \< 3) at investigated points from 
domain I. The numbers in parenthesis show the masses of squarks and gluinos (in GeV) 
at corresponding points. 3q 



